Geometric background charge: dislocations on capillary bridges. 
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Recent experiments have shown that colloidal crystals confined to weakly curved capillary bridges 
introduce groups of dislocations organized into 'pleats' as means to relieve the stress caused by the 
Gaussian curvature of the surface. We consider the onset of this curvature- screening mechanism, 
by examining the energetics of isolated dislocations and interstitials on capillary bridges with free 
boundaries. The boundary provides an essential contribution to the problem, akin to a background 
charge that "neutralizes" the unbalanced integrated curvature of the surface. This makes it favorable 
for topologically neutral dislocations and groups of dislocations - rather than topologically charged 
disclinations and scars - to relieve the stress caused by the unbalanced gaussian curvature of the 
surface. This effect applies to any crystal on a surface with non- vanishing integrated Gaussian 
curvature and stress- free boundary conditions. We corroborate the analytic results by numerically 
computing the energetics of a defected lattice of springs confined to surfaces with weak positive and 
negative curvature. 

PACS numbers: 45.70.-n, 61.43.Fs, 65.60.+a, 83.80.Fg 



When a hexagonal lattice is draped over a curved sub- 
strate, its bonds necessarily stretch or compress to ac- 
commodate the underlying curvature [TJ [2] . This phe- 
nomenon, known as geometric frustration, arises from 
the competition between local particle interactions, that 
favor an ordered lattice structure, and geometric con- 
straints that prevent local order from filling space [3UT2] . 
An intriguing connection has recently emerged between 
the frustrated order of twisted filament assemblies and 
the apparently distinct problem of curved crystalline or- 
der Q2]. 

A striking experimental realization of curved space 
crystallography is provided by colloidal monolayers con- 
fined on spherical droplets [TSJ [16] and on capillary 
bridges of varying positive or negative curvature [17] , see 
Fig. la. In both cases, the curvature-generated mechan- 
ical strains are screened by defects in the crystal lattice, 
especially dislocations or disclinations (see Fig. lb,c). 

In a spherical crystal, there must be an excess of twelve 
five-fold or (positive) disclinations simply because you 
cannot tile a sphere with hexagons, as demonstrated by 
a classic soccer ball. The presence of these twelve ex- 
cess disclinations is dictated by topology, energetics only 
fixes their positions at the vertexes of an icosahedron in- 
scribed in the sphere. The experimental and theoretical 
studies of Ref. [4j [15] have demonstrated that introduc- 
ing additional grain boundary "scars" emanating from 
each of the twelve positive disclinations becomes ener- 
getically favourable on large spheres, even if they are not 
required by topology. 

Defect nucleation on curved capillary bridges shaped as 
catenoids or barrels works differently than on the sphere: 
there are no topological requirements that force isolated 
disclinations in the ground state. Instead defects in the 
form of isolated dislocations are first nucleated to screen 
the curvature. Do "scars" exist on these surfaces even 
without the isolated disclinations from which they em- 
anate in spherical crystals? Recent experiments show 



clear evidence for the existence of finite length low-angle 
grain boundaries that do not originate from free disclina- 
tions [17] . These defect structures play an analogous role 
to pleats in fabrics - they accommodate for curvature by 
inserting and terminating extra rows of particles. 

The building block of 'pleats' is dislocations. In this ar- 
ticle we present an analytical and numerical study of the 
energetics of dislocations on capillary bridges of positive 
and negative Gaussian curvature. Our results highlight 
the ability of dislocations to screen curvature and pro- 
vide an accurate estimate of the forces acting on them 
as a result of the simultaneous presence of curvature 
and free boundaries. As we shall see, of crucial impor- 
tance is a subtle geometric force exerted on an isolated 
dislocation by a neutralizing "background charge" that 
counter-balances the effect of a non- vanishing integrated 
Gaussian curvature of the surface, thereby enforcing the 
stress- free boundary condition [5]. 

The elastic energy of a curved crystalline monolayer 
draped on a curved substrate can be written as (TJ [2] 

where Y is the Young modulus and x( x ) ls the Airy stress 
function that satisfies the biharmonic equation sourced 
by the Gaussian curvature and the disclinations 

y V\(x) = QaS(x - x a ) - G{x) , (2) 

a 

where q a and Xq, are the disclination charges and posi- 
tions respectively. 

The energetics of these elastic systems can be mapped 
onto a simpler electrostatic problem. Denote by x°( x ) 
the solution of Eq. ([2| in the absence of defects. The 
biharmonic equation can be solved in two steps [5j [9]. 
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FIG. 1: (a) Profiles of capillary bridges coloured by the mag- 
nitude of the Gaussian curvature on their surface (Eq.[9|. For 
the weakly curved capillary bridges shown here, the curvature 
G ~ R is approximately uniformly negative or positive. 
R z and are the radii of curvature in the vertical and az- 
imuthal direction respectively, (b) Disclinations in a hexag- 
onal lattice appear as five or seven-fold coordinated particles 
and carry an angular topological charge that can compensate 
the angular deficit created by the Gaussian curvature, (c) 
Dislocations in a hexagonal lattice appear as dipoles of discli- 
nations. The topological charge of disclinations appears in 
Eq. [2] as a source of geometric stress on the same footing 
as Gaussian curvature, suggesting a direct defect- curvature 
coupling, with five (seven) -fold disclinations coupled to posi- 
tive (negative) curvature. The energetics of disclinations are 
however more closely analogous to those of an electric charge 
that interacts with a geometric 'charge density' p, shown in 
panel (d). (d) p is the solution to a Laplacian equation (Eq. 
[3} with G as its source and an additional constant term pu 
that originates from the stress-free boundary condition. In 
the case of surfaces with unscreened curvature and free bound- 
aries, this term is akin to a background charge that neutralizes 
the unbalanced integrated curvature of the surface and leads 
to an energetic preference for the appearance of dislocations 
(dipoles) over disclinations, despite the uniform Gaussian cur- 
vature. 



First introduce an auxiliary function p(x) that satisfies: 

V 2 p(x) = G(x) . (3) 

The second step is to set p(x) as a source in the Poisson 
equation for x°( x ) 



1 

Y 



V 2 x G (x) = -(p(x)-p H (x)) , 



(4) 



where Ph(%) is a harmonic function that we will use to 
satisfy the boundary condition of vanishing stress at the 
edge of the capillary bridge [5j [9]. Note that the geo- 
metric charge density p(x) is analogous to a smeared out 
electrostatic charge that, as we shall see, interacts with 
the defects. 



In order to enforce the boundary condition, we need 
to determine the stress tensor o~ij(x) in an arbitrary set 
of curvilinear coordinates x — {x\,x<i\ which can be ex- 
pressed in terms of x( z ) an d the metric tensor g^ as [14]: 
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rg pr (d n d rX (z) - K r d qX {z)) 



(5) 



(6) 



The Christoffel symbols are denoted by T^. and 7^ = 
Cij/s/g is the covariant antisymmetric tensor, with the 
determinant of the metric g = ^11^22 • 

To facilitate comparison with our experimental studies, 
we focus on capillary bridges that are surfaces of revolu- 
tion whose radius r is approximately given as a function 
of the height z by: 



r(z) = + R z — R z cosh 



R z 



(7) 



This choice describes both the catenoids and barrels of 
Fig. la. Rc/y is always positive, while R z is positive on 
barrels and negative on catenoids. In both cases, the 
resulting metric reads: 



ds — cosh 



R z 



dz 2 



R(k + R z — R z cosh 



Rz 



(8) 

We assume that the total height of our catenoids and 
barrels is equal to H, so that — y < z < +y- The 
Gaussian curvature of these surfaces is given by: 

sech 3 (£) 



G(z) = 



R z (R<p - R z ) + R 2 Z cosh(^) 



(9) 



For the capillary bridges, we can calculate cr zz (z) and 
cr^z) from Eqs. ([Hp]) giving: 
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(R<f>-R z ) + Rz cosh (j^j 



dzX 
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cosh 2 ^) L 



1 z 

d z X~ ^ tanh (^) 9 *X 
rt z l\ z 



(10) 



(ii) 



The sign convention we adopt for is that the stress is 
negative if the material element is compressed and pos- 
itive if it is stretched. Intuitively we expect that when 
z ~ 0, near the neck of the catenoid (barrel), the applied 
stress is negative (positive). 

As we shall see, defect nucleation first occurs in the 
regime of small curvature H/R z <C 1, for which the cap- 
illary bridge can be viewed as a weakly deformed cylinder 
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FIG. 2: Energy (c,d) of a defected lattice of springs of unit rest length and unit spring constant, constrained to lie on the 
surface of capillary bridges (a,b). A dislocation, polarized so its five- fold particle lies above its seven- fold particle was inserted 
in the crystalline patch at positions of varying height in order to map out the potential felt by an isolated dislocation. The 
minimum energy of each configuration was found using a conjugate- gradient minimization similar to that used in ref [5] and is 
plotted (hollow squares in c(d)) for the capillary bridges shown in b(a). The flat-space energies (grey diamonds in c(d)) and 
energy of geometric frustration (black solid line in c(d)) of an un-defected patch are subtracted in (c,d) to reveal (red hollow 
circles) the predicted geometric potential (solid red line). 



of radius approximately equal to and height H. In 
this regime, we can Taylor expand all previous expres- 
sions in powers of z/R z to find a closed- form solution for 
X°( z ) m Eq. subject to the free boundary condition 
a? z (z = ±H/2) = 0. The result reads: 
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(12) 



where the last (constant) term, which is of no physical 
consequence, was chosen to set x G = on the boundary. 
Upon plugging Eq. ( [29] ) into Equations ( 10 ) and ( 11 ), we 
find (to leading order in H/R z ) that the stress cFf z ~ 
while <T00 « ?ar near the neck z ~ 0. These re- 
sults mean that there is spring compression for catenoids 
(R z < and cr9, < 0) and stretching for barrels (R z > 
and <T00 > 0). Note that the azimuthal component of 
the stress of geometric frustration, a9,(z) : changes sign 

for z > H/y/l2, a fact of crucial importance to deter- 
mine the geometric force exerted on dislocations by the 
curvature of the capillary bridge. 

According to standard elasticity theory, a dislocation 
in an external stress field, cr^(x), experiences a Peach- 



Koehler force, /(#), given by [18] 

fk(x) = e kj bi(Tij(x), 



(13) 



where b is the Burgers vector of the dislocation. Simi- 
larly, a dislocation introduced into the curved 2D crystal 
will experience a Peach-Koehler force as a result of the 
pre-existing stress field of geometric frustration cr^(x). If 

we choose b along its minimum-energy orientation for a 
catenoid (azimuthal, corresponding to a dislocation hav- 
ing its 7- fold defect closer to the neck) , we obtain a Peach- 
Koehler force f z (z) that points in the z direction and has 
magnitude given by: 



f z (z) = -ba%{z) 



(14) 



An explicit formula for f z (z) on capillary brid ges is ob- 
tained by plugging ~ d^x G ( z ) into Eq. (14) with 
the result 
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For a dislocation with its negative disclination facing the 
neck of a capillary bridge with negative curvature, the 
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force is repulsive (attractive), if it is located at a dis- 
tance z smaller (greater) than H/y/Vl. For a barrel, the 
prefactor R z switches sign, hence the minimum energy 
orientation for the dislocation is obtained for b < 0, eg. 
the positive disclination is facing the neck. 

The corresponding potential is obtained upon integrat- 
ing along the z direction (ignoring higher order metric 
effects) with the result 



minimum of <p(z) vs. G = 



Yb 



6R z Rrt 



H 2 



(16) 



Around the neck z > 0, the second term in Eq. (16) 
dominates and the potential <p(z) ~ 2^r H r z * s re P u l srve 
for the catenoid (Fig. Id) and attractive for the barrel 
(Fig. lc). For z < the dislocation acquires the opposite 
burger vector and the potential flips sign as expected 
given that the disclination closer to the neck has now the 
opposite sign. 

The second term in Eq. (16) is obtained from fixing 



the boundary conditions and can be traced all the way to 
the additional source pu in the Poisson equation (Eq. [4| . 
It is a contribution to the geometric force experienced 
by the dislocation that is linear in the Burgers vector 
and it is induced by the simultaneous presence of both 
the Gaussian curvature and the boundary. It should not 
be confused with the interaction of the dislocation with 
its own image which would be present also on a finite 
cylindrical patch. In the electrostatic analogy, the source 
pH subtracted from p(x) can be viewed as a neutralizing 
" background charge" which ensures that the area integral 
of the right hand side of Eq. Q vanishes. This condition 
forces the radial stress to vanish at the boundary. 

In order to test our analysis, we have performed con- 
strained energy minimizations of a hexagonal lattice of 
harmonic springs with the dislocation fixed at position z 
from the center of a patch of a catenoid (Fig. 2a) and 
a barrel (Fig. 2b). The gray hollow diamonds represent 
the energy of a patch containing a dislocation at different 
zona cylinder of the same size as the catenoid and bar- 
rel. It is very flat until the dislocation is placed close to 
the boundary and starts being attracted by its own im- 
age. This flat space contribution to the energetics purely 
depends on the presence of a boundary. We now subtract 
from the total energy (black hollow squares) the energy 
of the same patch in flat space (gray hollow diamonds) 
and the energy of geometric frustration of an un-defected 
patch (black solid line) and obtain (red hollow circles) 
for both catenoids and barrels the geometric potential, 
in agreement (red solid line) with that calculated in Eq. 
(16). The agreement with our analytical result (contin- 
uos line) that captures all curvature contributions to the 
energetics is very good. 

Figure 3b shows a systematic comparison of our nu- 
merical results (symbols) with the theoretical predictions 
(continuos lines) based on Eq. (16) for the z > portion 
of the geometric potential of a family of catenoids with 
increasing curvature. Figure 3a compares a plot of the 



R z Rd 



(the Gaussian curvature 



at the neck) determined numerically with the theoretical 
prediction. The agreement is very good, corroborating 
our analysis for a wider range of curvature than expected 
on the basis of our perturbative treatment that describes 
catenoids and barrels as small deformations from a cylin- 
drical geometry. The intuitive reason why the minimum 
of <p(z) is located at z ~ H/ \/l2, where changes sign, 
can be explained as follows. If the negative (positive) 
disclination is facing the neck of the catenoid (barrel), 
the extra row of atoms emanating from the dislocation 
is directed towards the z < z m i n (z > z m i n ) portion of 
the crystal that is stretched - this choice insures that the 
elastic energy is minimized. 

We can now use an electrostatic analogy to first deter- 
mine the geometric potential of a tight dipole composed 
of a positive and a negative disclination and then ver- 
ify that it indeed matches our geometric potential for a 
dislocation obtained analytically in Eq. ( [l6| ) and numeri- 
cally in Fig. 3. The energy of N disclinations labelled by 
the index a and with "topological charges" {q a = i 2 ^} 
reads j4l [9] 
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(17) 



where is the Green's function of the biharmonic op- 
erator. The source Mix) is given by 



M(x) 
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q a 8(x, x a ) — G(x) . 



(18) 



Equation (17) implies that the potential U(z), experi- 
enced by an isolated disclination on a curved substrate 
with Gaussian curvature G(z) satisfies the following bi- 
harmonic equation 



V A U(z) = -G(z) 



(19) 



and is therefore proportional to the Airy stress function 
X G (z). 

The geometric potential of a dislocation can be con- 
structed by calculating the energy of a disclination dipole 
whose moment qdi = e^bj is a lattice vector perpendic- 
ular to b that connects the two points of 5 and 7-fold 
symmetry along the surface. Upon taking the product of 
d times the gradient in the z direction of the disclination 
energy qU(z), including the factor v /g^ = cosh(z/r 2 ), 
we obtain the potential <j>{z) and upon taking a second 
gradient the geometric force f z (z) 



fz{z) = 



1 



qd 



d z U{z) 



(20) 



Once the identification U(z) = X ( z ) an d qd = b is 
made, Eq. (20) becomes equivalent to Eq. ( [l4| ) . The 
energy qU(z )of a disclination with topological charge 
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FIG. 3: Geometric potential for an isolated dislocation on 
surfaces of increasing Gaussian curvature, obtained in a sim- 
ilar way as the curves shown in Fig. 2. The scaling of the 
minimum energy as a function of the Gaussian curvature at 



the neck G — R R ^ , is shown in (a), while a comparison 

of predicted geometric potentials (dashed lines) and energies 
found numerically (solid rhombi) is shown in (b) 



q is readily obtained from solving the Poisson equation 
V 2 U(z) = - (p(z) - ph(z)), see Eq. Positive (nega- 
tive) disclinations are repelled (attracted) from the neck 
of the catenoid (barrel) at z ~ by the integrated back- 
ground source p(z) — Ph(z). 

We can also determine the geometric potential X(z) of 
an interstitial as the product of the local area expansion, 
5 A, generated by the interstitial, times the local pressure 
p(x) = — V 2 x G ( x ) generated by the curvature. The pres- 
sure (or equivalently the interstitial geometric potential) 
is proportional to the right hand-side of Eq. Q, eg. it 
is equal to the source of the Poisson equation including 
the crucial boundary term pu{x). To leading order in 
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FIG. 4: Numerically calculated energy (black hollow circles) 
and predicted energy (dashed line) of an interstitial particle 
on a capillary bridge having negative (a) and positive (b) 
curvature. The energy of an interstitial is given by p(z)b 2 and 
therefore provides a means to directly probe the geometric 
potential p(z) on a surface. 



the surface deformation we find 
Yb 2 
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(21) 



where we assumed 8 A ~ b 2 . 

The analytical expression for X[z) derived in Eq. (21) 
is plotted as a continuos line and compared to numerical 
results for catenoids and barrels in Fig. [4j Note that the 
only effect of switching the sign of the Gaussian curva- 
ture in going from catenoids to barrels is the relative sign 
of R z and R^. This determines the overall negative (pos- 
itive) pre-factor of the parabolic potential for catenoids 
(barrels) in Eq. (pTl) . 



The numerical data shown in Figure [4] are obtained 
from constrained energy minimization algorithms simi- 
lar to the ones used to determine the dislocation poten- 
tial. As shown before, we first determine numerically 
the interaction of the interstitial with its own image on 
a fiat cylinder, and then subtract it from the total en- 
ergy to isolate X[z) for a catenoid (panel a) and a barrel 
(panel b). The agreement is very good confirming the 
intuitive expectation that interstitials are repelled (at- 
tracted) from the neck of catenoids (barrels) where the 
frustrated lattice is compressed (stretched). Exactly the 
opposite trend is expected for vacancies which represent 
an area deficit. 

As we have illustrated, charge neutral dislocations can 
screen the Gaussian curvature of the capillary bridges as 
effectively as disclinations. We can estimate the thresh- 
old integrated Gaussian curvature required to make dis- 
locations nucleation energetically favorable. On capillary 
bridges, this instability is a precursor to pleating: the for- 
mation of charge neutral grain boundary scars observed 
in the experiments of Ref. [17 . In what follows, we 
capture the onset of this instability by calculating the 
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energetics of a pair of dislocations with opposite Burg- 
ers vectors oriented so that the 7s face the neck of the 
catenoid. 



Upon using Eq. (14), the difference in energy, AE(z), 
between the curved crystal with the two dislocations at 
positions ±z and the curved crystal without any defects 
is readily estimated 



AE(z) 
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2b 



b (z')dz' . (22) 



The first term is the flat-space logarithmic interaction 
energy between two dislocations being pulled apart from 
the neck which approximates the interaction between dis- 
locations to zero order. The second term accounts for the 
interaction of each of the two dislocation with the curva- 
ture. The energetic condition for the nucleation of a pair 
of oppositely oriented dislocations is AE(z) < 0. 

Upon plugging into Eq. (22) we obtain for AE « 

0: 



H 2 



9v^3 RqRz 



b (H 



(23) 



The Burgers vector b is typically of the order of the lattice 
spacing a leading after simple rearrangements to 
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(24) 



The energetic condition can now be rewritten as a ge- 
ometric condition on the threshold integrated Gaussian 
curvature: 
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(25) 



Note that for large samples the threshold integrated cur- 
vature is lowered: the instability is not suppressed but 
rather enhanced. 

The geometric criterion for the nucleation of a disloca- 
tion pair, derived in Eq. (25), is qualitatively similar to 



heuristic estimates of the onset of pleating. Both insta- 
bilities are triggered by the energetic advantage derived 
from adding (or subtracting) rows of atoms where the 
curved crystal is most stretched (or compressed). Nucle- 
ating both dislocation pairs or pleats is associated with 
an energetic penalty that diverges slowly as the loga- 
rithm of the system size. By contrast, the energy Ef of a 
stretched or compressed curved crystal characterized by 
a bond angle rotation 6 reads 

E f ~ Y0 2 H 2 . (26) 

Note that it diverges quadratically with the size of the 



system. 

A single pleat is essentially a low angle grain boundary 
of finite length H whose energy can be estimated as [18] 



E p ~ -YbH 0\n0 



(27) 



As long as the curvature induced angular stretching (or 
compression) is small, it is energetically advantegeous 
to stretch (compress) bonds in the crystal rather than 
adding extra rows of atoms because Ef ~ 2 < E p . 

Pleating becomes energetically favorable above a crit- 
ical 6 C 



H 



In Or. 



(28) 



Equation (28) can be solved self-consistently. Neglecting 



doubly logarithmic corrections leads to 



b (H 



(29) 



The curvature induced angular stretching (compression) 
is approximately equal to the integrated Gaussian curva- 
ture. This gives the criterion 



GdA 



I'M- 

H \ a 



(30) 



where the Burger vector is approximately one lattice 
spacing long. Note that Eq. (30) provides an estimate of 
the same order as Eq. (25). 

In this article, we have demonstrated that continuum 
elasticity can capture the energetics of defects on curved 
capillary bridges with free boundaries, as illustrated by 
the favorable comparison between our analytical results 
and numerical energy minimizations of discrete spring 
models. Crucial to this agreement is the inclusion of a 
subtle geometric interaction between defects and a neu- 
tralizing background charge generated by the boundary. 
Our results are in agreement with the experimental find- 
ings of Ref [17] and show that charge neutral pleats or 
dislocation pairs are effective in screening the Gaussian 
curvature on curved capillary bridges. 
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